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ABSTRACT 


A technique of statistically adjusting dynamical fore- 
casts of tropical cyclone motion was tested. All tests were 
performed with operationally-analyzed data from the U.S. Navy 
Fleet Numerical Weather Central (FNWC). Three sets of regres- 
Sion equations were developed to modify forecasts of typhoon 
maaGkS. The first=set of equations was based only on forward 
intejration of the FNWC Tropical Cyclone Model (TCM) for 28 
cases in 1975-76. An independent sample of cases from 1977- 
78 indicated that the first equation set was based on too 
Small a sample size, especially considering the anomalous 
nature of the 1975-75 storm tracks. A second equation set 
based only on forward integration of the TCM was derived from 
ees COunmuacemno@acdsts from 1975-78. Results from the ex- 
periments with these equations indicate that systematic data 
Mice mOnhemecmuamomean be wsed to statistically adjust forecast 
Storm tracks. The second equation set based on forward inte- 
gration showed improvement over the unmodified model predic- 
tions at all forecast times. A third equation set based on 
forward and backward integration of the TCM exnlained the 
greatest amount of variance of all the equation sets. Ina 
dependent test of these equations using 31 of the 1977-78 
Cases olemeemmiavvmy th Fleet error goal of 100 and 150 
nautical miles at 48 and 72 hours was nearly met. 
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I. ANGER OME TT GN 


The tropical storm is one of nature's most destructive 
phenomena, making typhoon prediction a major concern world- 
wide. The need to forecast storm conditions has produced a 
Variety of prediction methods. Subjective methods, includ- 
ing persistence, have generally produced fairly good short- 
range forecasts. An especially difficult problem for any 
Subjective or objective prognostic scheme is the occurrence 
of storm recurvature. This problem was significant in the 
western North Pacific Ocean during the 1975 typhoon season, 
and to a lesser extent during the 1976 typhoon season. An 
unusual number of storms tracked northward, while others re- 
curved, As a result, the forecast error statistics were 
higher than normal (Annual Typhoon Report, 1975). The 
United States Seventh Fleet Commander, noting the necessity 
to improve the forecast errors, has levied a requirement for 
tae Joint Typhoon Warning Center (JTWC) to achieve maximum 
maimecasteerrors of 50, 100 and 150 nautical miles for 24, 
momande,e Wwours, respectively. In comparison, the 19/7 JTWC 
average position forecast errors for tropical cyclones at 24- 
48-, and 72-hours were 140, 266 and 390 nautical miles, re- 
spectively (Annual Typhoon Report, 1977). The long-term 
error trend, as indicated by the five-year mean error, has 
been increasing since 1972. In 19/77 the operational [CM 
mOduced mean vector errors of 138, 262 and 450 nautical 


miles for all tropical cyclones (Annual Typhoon Report, 1977). 
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Although progress has been made in the numerical simula- 
faceoff ChOpical Storm characteristics such as intensity and 
spiral bands (Anthes et al, 1974), the major emphasis has 
focused on storm motion (Hovermale et al, 1976; Ley, 1975). 
Multi-level, nested-grid models are being developed in an 
attempt to improve forecast positions over those of analog, 
persistence and statistical techniques. In spite of the 
sophistication of some numerical models, undesirable conse- 
quences may result from initialization with poor or limited 
data. One test of the dynamical models is the accuracy of 
mcm nbcudiestorm track. If the short-term forecasts are 
inaccurate, a dynamic model cannot be expected to produce 
accurate, extended forecasts. Some dynamic models (Hovermale 
edemeooe Ley, 1975; Hodur and Burk, 1978) predict tracks 
with systematic bias relative to the actual track and com- 
fomiey predict motion which is too slow. A likely source 
of error seems to be in specification of the initial data 
fields, although inadequate resolution in the numerical 
moe! could also explain the biases. 

Since the resources required to increase high quality 
data or to run sophisticated numerical models are expensive 
and time consuming, alternative methods for improving model 
forecasts should be considered. Statistical-dynamical 
schemes for predicting tropical storm motion (Renard et al, 
Mp2 erely heavily on current storm motion... If position 
GENOrtsMOL tropical storms are not timely, then the statisti- 
cal schemes may be impaired. An approach which circumvents 


many of these problems is to use the model forecast positions 
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themselves to statistically modify the predicted track posi- 
mons lt 15 the primary objéctive of this thesis to de- 
velop and evaluate statistical regression equations for 
adjusting dynamically predicted storm tracks from the Tropi- 
cal Cyclone Model (TCM) used by the U.S. Navy Fleet Numerical 
Weather Central (FNWC). 

The basic model used for these experiments is the primi- 
tive-equation, three-layer, tropical cyclone model developed 
by Elsberry and Harrison (1971) and Harrison (1973). Although 
this model is capable of triply-nested operation (Harrison, 
Mewes, results from Ley and Elsberry (1976) show that the 
coarse and nested grids produced nearly identical results in 
a selected case study based on hand-analyzed data. In 1975 
FNWC adapted a coarse-mesh (Ne three-layer, dry version of 
this model for tropical cyclone prediction. Preliminary 
results with operational data were presented by Hinsman 
mere tor 1975 and 1976 data. For the 1977 typhoon season, 
the model was modified to include a biasing technique sug- 
gested by Shewchuk and Elsberry (1978). The use of the 
forecast stream function field to determine the cyclone 
forecast positions reduced errors in the relative vorticity 
tracking used prior to this modification (Shewchuk, 1977). 
The current operational version of the model has boundary 
conditions which are insulated, free-slip walls on the 
morth and south and cyclic on the east and west. Since it 
was Suspected that these boundary conditions could adversely 
Geer tie forecast storm track, Hodur and Burk (1978) in- 


corporated one-way interactive boundaries. Substantial 


ps 








improvement in storm track prediction was made compared to 
the previous channel version of the model, The open 
boundary version of the TCM was used in the experiments 


described herein. 
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IIT. BACKWARD INTEGRATION 


fee ide MODEL 

A detailed description of the Tropical Cyclone Model 
(TCM) is presented in Appendix A. The TCM is a three-layer, 
primitive-equation model in pressure coordinates. It is 
an open boundary model with one-way interactive boundary 
Sinciinonseone tie worth and South and cyclic conditions on 


the east and west boundaries. 


B. MODEL MODIFICATIONS 

ailieenorward wntegration of the model was carried out 
uSing the version of the TCM described by Hodur and Burk 
(1978). The backward integration was carried out using a 
negative time step of -600 seconds unless the northern 
boundary location was greater than 40°N. In the latter 
case, the negative time step was reduced to -450 seconds 
after the ori NovedwassUtirie€ilent distance to the north. 

In the forward integration, heat was added to the storm 
center as defined by a minimum wind at 1000 mb. The opur- 
Mose OTetnenneating function 18 to counteract the dispersion 
of the vortex due to the finite differencing (Ley and 
Elsberry, 1976). If heat had been added during the back- 
ward integration, this would have contributed to a better 
definition of the storm center. However, this presents a 
Physically unrealistic situation for a typhoon moving back- 


ward in time. For this reason, the heating function was 


e/ 









Set equal to zero for all backward integrations. Negative 
heating was rejected because it led to premature dispersion 


Sieene storm, which became difficult to track. 
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ITI, APPROACH TO TRACK MODIFICATION 


Two types of regression equations were tested here. One 
type will have predictors based only on forward integration 
of the TCM, while the other set of predictors will be based 
On both forward and backward integration of the TCM. The 
open boundary version of the TCM (Hodur and Burk, 1977; Hodur 
and Burk, 1978) was run with 46 cases (9 storms) from 1975 
and 1976. All cases were based on oOperationally-analyzed 
data obtained from FNWC. The regression equations developed 
from this sample were then tested against independent cases 
in 1977 and 1978. A second setof regression equations based 
only on forward integration of the TCM are derived using the 
combined data sets of 19/5-78. 

The final set of regression equations are based on forward 
integration of the model to 72 hours, as well as backward 
integration for 36 hours. The backward integration should 
reveal the effects of systematic model and data errors. 
Fundamental assumptions are that the model-related errors 
menGg to De systematic, and adjustment for data errors is 
possible where, in the absence of observations, the initial 
analysis reverts to the east-west flow appropriate to the 
climatology. Backward integration will increase the number 
Of predictors available to explain the variance between ICM 
forecasts and best track positions. This should lead to 
improved regression equations witn a higher explained variance. 


The predictors using backward-integrated positions can be 





used to adjust the forward-integrated positions using the 
Same initial data. 

Both the forward- and backward-predicted tracks were 
compared to the best tracks at their corresponding times. 
Next, regression equations to adjust the predicted track to 
the actual track uSing stepwise regression were generated 
fee tne “Statistical Package for te Social St@iences (SPSS). 
The regression equations were developed to make corrections 
acyl cenours up to /2 hours. If a TCM run was incomplete, 
alternate sets of equations were derived according to the 
length of the TCM run. The options considered were 36-, 48-, 
60-, and 72-hour TCM runs, with and without backward integra- 
Eon. 

If successful, the advantage of this approach would be 
that use of simple regression equations would require much 
less computer time than more sophisticated dynamic models. 
mus, it may be possible to produce tracks that are more 
accurate than ‘the ordinary open boundary TCM without a 
large increase in computer resources. An operational advan- 
tage of this approach is that no warning positions are 
required for computation of any predictors in the regression 
equations. The method and results of this approach to storm 


track prediction are presented in the following sections. 
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IV. DEVELOPMENT OF THE REGRESSION EQUATIONS 





A. THE PREDICTORS AND PREDICTANDS 
Using TCM forecast runs versus best tracks (Annual Typhoon 


Reports, 1975-77), 12 predictands were derived by computing 





the east-west and north-south differences between positions at 
corresponding times (see Fig. 1). Storm positions at each 
forecast time are adjusted by two regression equations, one 
for the east-west direction and the other for the north-south 
feeection, hhus, for a /2-hour TCM forecast run, a total of 
12 regression equations would be used to modify storm tracks 
iii 2-nour increments from 12 to 72 hours. 

Predictors used in these equations were model-predicted 
displacement and velocity, broken into components along the 
east-west and north-south directions. The Julian day, lati- 
made and iOngitude of the initial position of each TCM fore- 
cast run are also included as predictors in the regression 
equations. A schematic illustration of the intervals over 
which the predictors were calculated for the forward integra- 
mon runs 1s Shown in Fig. 2. A complete 11St om predictands 
and predictors, with the times for which they were computed, 


appears in fable l. 


Eee merwOD OF EQUATION DERIVATION 
The regression equations were derived using the Statisti- 
Egimeedekage for the Social Sciences (itie, et al, 1975). 


Cases with missing values for predictors or predictands were 
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Figure 1. Depiction of the model errors (i.e., the predic- 
tands), which are the difference between the best 
tyack and forecast positions, which are shomn 
above as AXi o> AY 9. AX oq» AY 54> etc. 
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Figure 2. Depiction of intervals over which displacements 
and average speeds were computed using only for- 
ward integration of the open boundary TCM. Un- 
modified TCM forecast positions (0) are shown at 
12 hour intervals. 
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TABLE I 





Predictors/predictands used to develop regression equations 
| for typhoon track modification based only on forward inte- 
| gration of the open boundary TCM. 

imeeePredictands: AX, AY 


Times at which predictands are computed: 


ieee. 2O, Same O00 ge 2 hins 





wee rredivctors: ig A) (AURA) 
Time intervals over which each predictor was computed: 
00-12, 12-24, 24-36, 36-48, 48-60, 60-72, 
00-24, 12-36, 24-48, 36-60, 48-72, 00-36, 
00-48, 00-60, 00-72 hrs 


feeeiti tia! Position Predictors: Julian Day, Latitude, 
Longitude 


Times: 


Initialization time of a given TCM forecast run 
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automatically eliminated from all calculations. Such cases 
arise because the tracking routine is not always able to 
follow the storm center throughout the 72-hour interval. In 
the forward integration tests of Hodur and Burk (1978) with 
the open-boundary conditions, only 28 of the 46 cases extended 
moe? 2 hours. 

To avoid the problem of uncomputable predictors when 
the duration of TCM runs was less than 72 hours, it was 
decided to derive alternate sets of equations based on the 
duration of the TCM forecast run. The forecast lengths con- 
Sidered were 36, 48, 60 and /2 hours. For example, if the 
TCM produced only a 48-hour storm track forecast, no regres- 
Sion adjusted track would extend behond 48 hours, and only 
predictors in the range 00-48 hours were considered when 
deriving the regression equations to be applied to a 48-hour 
fen run. 

All predictors listed in Table I were considered for in- 
clusion in each regression equation. Selection of predictors 
Was stopped when the next variable in the stepwise regression 


explained less than 1% of the variance. 


Pee FIRST EQUATION SET BASED ON FORWARD INTEGRATION 
fomemlo75-76 cases, as well as those in 1977 and 19/78 to 

be discussed later, are based on cases with relatively well- 

developed storms which are more readily modeled by the coarse- 

mesh TCM. It should be noted that the operational version of 

the TCM is used only for storms exceeding 50 knots. Inclusion 


of weaker tropical storms would likely increase the variance 


Zo 








between the forecast and best track positions, thus making 
it more difficult to derive stable regression equations, f 
sample equation for the regression adjustment along the x- 
axis at /2 hours is shown below: 
DXER72=-537.0151+39.2088(XXLAT)+68.1605(VX1224) 


- 50.1251(VX0024)-43.8154(VY0048) 
- 1.4956(DX6072)+ 0.6102(JULDAY) 








Velocity was the most Frequently -selected preadpetor. [hs 

| Mesuit 1S not too surprising if one recalls that a common 
mort eof dynamic models is prediction of motion which is too 
Slow (Hovermale et al, 1976; Ley and Elsberry, 1976). The 
version of the TCM (Hodur and Burk, 1978) used in these ex- 
Semen ts 1S nO exception. Predictors with their associated 
time intervals beginning, ending or overlapping the valid 
time of a given regression equation were often selected. In 
General, this suggests a sort of "Statistical extrapolation” 
in which past, present and future forecast motions are used 
to correct the forecast track. 

The average explained variance of the regression equations 
appears in Table II. This parameter measures the strength 
of the linear relationship between the multi-linear regres- 
Sion equation value and the observed value. The amount of 
explained variance in Table II generally decreases as the 
TCM runs become shorter in duration and the number of avail- 
able predictors is thereby reduced. In general, longer TCM 
forecast runs permitted use of regression equations which 


corrected for more of the variance from the best track. 
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TABLE 11 


Average explained variance of the regression equations. 


Forward Integration (1975-76 Cases) 


feoene TCM Run 81.6% 
60 Hr TCM Run 77.6% 
48 Hr TCM Run 75.4% 
36 Hr TCM Run 68.4% 





Forward Integration (1975-78 Cases)’ 


| teense TCM Run 2) See: 





) Aomadr TCM Run 16.5% 
48 Hr TCM Run A bee 
B® 636 Hr TCM Run 43.2% 


Forward and Backward Integration (1977-78 Cases) 


[oor TCM Run 86.2% 
Boene TCM Run 84.3% 
48 Hr TCM Run 83.8% 
po Hr TCM Run 82.3% 
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Shortcomings (see Fig. 3) of the one-way (OW) interac- 
tive boundary Version of the TCM are that the forecast 
tnack typically runs to the left of the storm path and the 
velocities are generally too slow (Hodur and Kurk, 1978). 
These systematic errors of the TCM seem to contribute 
Significantly to the large amount of explained variance in 
the regression-produced typhoon positions (see Table I!). 

It should be recalled that the 1975 typhoon season, 
and to a lesser extent the 19/76 season, experienced a high 
frequency of storm recurvature as well as a large number of 
storms which tracked northward. This bias was present in 
the data sample used to derive the regression equations and 
had considerable impact on the selection of predictors and 


the computation of constants and coefficients. 


D. SECOND EQUATION SET BASED ON FORWARD INTEGRATION 

To increase the number of typhoon cases used to derive 
the regression equations, TCM forecast runs from 1977 and 
1978 were included in the data set. These ICM forecasts 
were originally used to make an independent test of the 19/75- 
76 regression equations. However, the majority of the 
typhoons in 1975 and many in 1976 had significant recurvature. 
The inclusion of the 197/-78 storm tracks should make the 
Sample more representative of the various tracks found in 
tne western North Pacific Qcean area, 

Because some TCM runs did not extend to 72 hours, only 
61 of 90 cases were available for derivation of the regres- 


Sion equations. As with the first set of equations based 
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Figure 3. Forecast tracks of Typhoon Marie produced by the 
OW model compared to the 8-13 April 1976 best 


track poSitions. “Baen cCiretle represents a 12ahiour 
increment (after Hodur and Burk, 1978). 
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on the forward integration of the TCM, alternate sets of 
equations were derived based on 72-hour, 60-hour, 48-hour 
and 36-hour TCM runs. When less than 1% of the variance was 
explained by a variable, selection of predictors was halted. 
The complete set of regression equations is listed in 
Appendix B. 

The expanded data sample (61 cases) contained a greater 
variety of storm tracks and reduced the recurvature bias of 
the 1975-76 typhoon cases. Storm tracks in the 1975-78 
data set used in these experiments are characterized by four 
general categories: westward, northwestward, northward, 
and recurving paths. With this greater variety of storm 
paths, the amount of explained variance of the regression 
equations (see Table II) dropped commensurately. For example, 
the regression equations for a /72-hour TCM run incurred 
approximately a 30 percent decrease in explained variance. 
With the exception of a small variation at 48 hours in Table 
Pi, tie general trend was again a reduction of explained 
variance as the duration of the TCM forecast runs decreased 
meeoomnours. AS discussed previously, this characteristic 
is attributed to having fewer predictors available to explain 
the variance when the TCM forecast is of shorter duration. 
These results reaffirm that this method of adjusting TCM 
forecasts is at its best when the TCM forecast extends to 
Famnaurs, rather than a shorter forecast interval. 

Velocity was again the most frequently selected type of 
predictor, indicating an attempt to compensate for the over- 


all slowness of the TCM. However, there was a better 
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balance of velocity and displacement predictors in the 1975- 


78 equations as compared to the 1975-76 equations. 


E. EQUATION SET BASED ON FORWARD AND BACKWARD INTEGRATIONS 

The third set of regression equations was based on storm 
positions derived from both backward and forward integration 
of the TCM. Additional velocity and displacement predictors 
based on backwards integration of the TCM as indicated in 
Fig. 4 and Table III were added to the data set. These pre- 
dictors are computed in the Same manner as those based on 
forward integration. All predictors listed in Table I and 
Table III were considered for inclusion in each regression 
equation. Regression equations were again developed for TCM 
momecast runs of 36, 48, 60 and 72 hours duration. 

The regression equations which included backward-inte- 
grated positions had the highest average explained variance 
(see Table II). As TCM forecast runs become shorter in 
duration, the amount of explained variance in the equations 
decreases due to fewer predictors being available (see 
Table II). The complete set of regression equations is 
listed in Appendix B. 

The regression equations took advantage of the systematic 
errors inherent in the initial fields and in the numerical 
model. The most favored predictors from the backward inte- 
waauione occurred in the interval from -12 to 00 hours. For 
the 12 equations used when the TCM was integrated forward 
mom eehours and backward to -36 hours, predictors in the 


interval -12 to 00 hours appeared in 9 equations. The 


3H 





7276) 


60 O 


48 O 


Figure 4. Depiction of intervals over which displacements 
and average speeds were computed using bo Gi 7 0. - 
ward and backward integration of the open bound- 
ary TCM. Unmodified TCM forecast positions (0) 
are shown at 12-hour intervals. 
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TABbe lil 





| Additional predictors used to develop regression equations 
| for typhoon track modification based on the backward inte- 
gration of the open boundary TCM. 


Preaicetors: Ax, AY, u, v 
Times = O0O-M12, Mi12-M24, M24-M36, 12-M12, 


OO-M24, M12-M36, OO-M36, 24-M24 hrs 


where 00 = initial time 
M24 = minus 24 hours from initial time 
12 = plus 12 hours from initial time 
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velocity along the y-axis from -12 to 00 hours appeared in 7 
of the equations, andin each case it explained the most 
Variance of any aed Leer in the equation. Approximately 
35 percent of all the predictors in these 12 regression 
equations were computed from backward-integrated storm 
positions. Velocity was the most frequently selected pre- 
dictor from all forecast intervals. As will be indicated 
in later examples, the Fegresé1 on se@uaimiomne appeared to be 
compensating for the generally slow motion of the TCM. It 
should be noted here that these regression equations and 
the sample cases which follow were based on a limited data 


Sereotf ol storm cases from 1977 and 19/78. 
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Ve RESULTS 


A. FIRST EQUATION SET BASED ON FORWARD INTEGRATION 

The primary purpose of the tests with the 1975-76 cases 
was to make a preliminary evaluation of the regression 
equations which were derived uSing this data set, Because 
this data set was used to derive the regression equations, 
these experiments are referred to as a dependent test. 

A sample of adjusted storm positions based on the first 
equation set is shown in Figs. 5-7. Results with Typhoon 
June shown in Fig. 5 and Typhoon Marie depicted in Fig. 6 
were very encouraging, Typhoon Marie is the same storm 
which appears in Fig. 4. Storm recurvature was accurately 
depicted, and storm velocity was markedly improved in these 
two storms. Note, however, that the velocities in the 
regression tracks are generally greater than the best track 
meocwties. This characteristic of the regression equations 
appears to be an attempt to compensate for the slowness of 
fiemncM. In the majority of cases in this sample, the 
storms were excessively accelerated due to this feature of 
the equations. Since the regression equations were based 
Primarily on predictors derived from TCM forecasts, it is 
clear that the goodness of the regression-adjusted storm 
track is dependent on the quality of the TCM forecast itself. 
The adjusted track for Typhoon Rita (August 1975) shown in 
Fig. 7 illustrates how a poor TCM forecast will lead to an 


extremely radical regression adjustment. This example 
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maaicates how small TCM velocities, as well as radical 
changes in the forecast storm track, will result in a 
regression-adjusted track with extremely high velocities 

and unrealistic variations in the storm path. This behavior 
of the regression equations seems to indicate that obviously 
erroneous regression positions can be used as a basis to 
reject TCM forecast positions as well. 

The mean forecast errors of the 1975-76 dependent test 
Sample are listed in Table IV. The sample improved on the 
forecast errors of both JIWC and the unmodified TCM fore- 
@asts for, 19/5-76. At 72 hours with a sample size of 28, the 
regression equation errors were also less than the U.S. 
Navy 7th Fleet goal of 150 nm. However, it is not expected 
that this pattern would be repeated in an independent sample 


of typhoon cases. 


Pee NNDERENDENT TEST OF FIRST EQUATION SET 

An independent sample of 44 cases from the 1977-78 typhoon 
seasons was then used to test the regression equations de- 
rived from the 1975-76 cases discussed in the previous 
SeGulanemers indicated by the statistics in Table V, the 
results of the independent test were generally very poor. 
In each category the modified TCM tracks were worse than the 
unmodified tracks. The large errors are attributed to the 
unstable nature of the regression equations from the small 
and very homogeneous sample of anomalous storm tracks used to 
derive the regression equations. This led to excessive 


velocities and erratic tracks in the regression-adjusted 
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Secu positions for 1977-78. When thie TOM forecasts had 

large deviations from the best track, the regression positions 
incurred commensurately larger errors. [It was clear that 

the number of TCM forecast cases had to be increased as much 
as possible to increase the stability of the regression 
equations. 

The regression equations used in these tests did not 
improve storm positions by taking advantage of the systematic 
errors of the TCM, which is often too slow and to the left 
of the best track. Typhoon Lola is depicted in Fig. 8 and 
mmanredtes Slight corrective shifting of the regression posi- 
moms to the right of the TCM forecast track. HoWever, over- 
compensation iS again evident in the velocity components. 
More typical of the independent test results is the adjusted 
meackenror Typhoon Lucy shown in Fig. 9. Velocity errors were 
large and the track was erratic. In many cases in the inde- 
pendent test, an erratic track, such as the regression- 
adjusted path for Typhoon Lucy, occurred when the unmodified 
TCM forecast track also incurred large errors. The single 
redeeming feature of these 19/75-/6 equations may be that 
highly erratic regression positions could be indicative of a 
poor TCM forecast. This information may be useful in leading 
the typhoon forecaster to reject both the modified and the 


unmodified TCM guidance in making his decision. 


C. SECOND EQUATION SET BASED ON FORWARD INTEGRATION 
The second set of regression equations was obtained from 


the combined storm cases in the 1975-76 and 19/7/-78 samples. 
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Because this exhausts the sample, these equations can only 
be evaluated as a dependent test. 

This second os of equations, with a larger sample of 
61, was more stable overall than the first set of regression 
equations. While the regression-adjusted positions still 
Increase the velocity excessively, the apparent result is 
that this overcompensation is decreased in the 1975-78 re- 
gression equations versus those for 1975-76, 

The five typhoon cases shown previously in Fig. 5 through 
moe appear again in Fig. 10 through Fig. 14, but this time 
illustrating the regression-adjusted positions of the 1975- 
78 equations. The track for Typhoon June shown in Fig. 10 
was not adjusted as well in this test, but it serves to 
illustrate that the regression adjustments cannot always com- 
pensate for the tendency of the TCM to not predict recurvature. 
Presumably if the TCM track had been more northwesterly, the 
adjustment would have been toward more recurvature, The 
track for Typhoon Marie depicted in Fig. 11 still indicates 
a reasonable adjustment for recurvature, but has incurred a 
marked decrease in velocity adjustment. Although the adjusted 
speed of movement for Typhoon Rita shown in Fig. 12 is slower 
mianeihneat on Fig. 7, the track 1S not significantly changed 
from the previous result. However, this case again suggests 
that a radical and obviously erroneous regression adjustment 
May serve as a basis to reject the TCM forecast also. Typhoon 
Lola (September 1978), illustrated in Fig. 13, has a much 
improved regression adjustment of the storm track compared 


with Fig. 8. The velocity corrections in this case are much 
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closer to the actual velocities which appear in the best 
medek positions. The adjusted track for Typhoon Lucy is 
saownm In Fig. Ill nie 1s to be compared wrth Fig. 9. The 
regression-adjusted positions are still poor, but not 
nearly as erratic. Note that the regression-modified posi- 
tions do suggest some recurvature, although inaccurately, 
and that the decreased velocities are much more realistic 
with the 1975-78 regression equations, 

Generally, the results of these experiments suggest 
that a further increase in the number of TCM forecasts in 
the data sample would lead to a further improvement in the 
regression-modified storm positions. Qverall, the errors 
incurred by the regression-adjusted tracks were less than 
those of the unmodified TCM storm path (see Table VI). The 
improvements are especially noteworthy at 48 hours and 72 
menus - 

In forecasting typhoons which tracked westward, the un- 
modified TCM forecast positions usually proved to be very 
accurate, with the regression-modified path seldom improving 
eoetne TCM, One case (Typhoon Rita, October 1978) in which 
the regression equations improved the forecast is shown in 
pice lS, the improvement 1S the result of an increase in 
storm velocity by the regression equations. Several other 
examples of enhanced track predictions are shown in Fig. 16 
through Fig. 18. These examples show some improvement in 
direction, but the improvement in speed of movement is 


especially rewarding. 
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D. EQUATION SET BASED ON FORWARD AND BACKWARD INTEGRATION 

These experiments with the 1977-78 typhoon cases were 
performed to test the regression equations derived from the 
1977-78 data set which included backward-integrated posi- 
moms tO =36 hours. As the complete set of backward tracks 
was used in the derivation of the regression equations, 
these experiments refer only to the dependent sample. 

The results of the dependent wane with forward- and 
backward-integrated positions were the best of all experi- 
ments conducted during this research. This equation set in- 
curred the lowest mean errors overall (see Fig. 19). The 
regression-modified storm position errors are very close 
to meeting the U.S. Navy 7th Fleet error goal (see Table VII). 

Weacks for Typhoons Lola, Rita, and Babe shown previowsly 
mumeree. «1S, 15, and 1/7, respectively, appear again tn Fig. 20 
through Fig. 22 with adjusted positions based on forward- 
backward integration. Typhoon Lola (September 1978), shown 
in Fig. 20, is more erratic than the same case (Fig. 13) 
using regression equations based only on forward integration, 
However, the regression positions that are based on forward- 
backward integration are fluctuating on either side of the 
best track. The variation about the best track is attributed 
primarily to the small size of the data set (31 cases) used 
to derive the equations. This "Saw-tooth" variation about 
the best track is characteristic of the behavior storms which 
tracked westward in this experiment. As in previous tests, 
it proved difficult to improve on the unmodified TCM forecast 


positions when typhoons were tracking approximately westward. 
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Figure 19. Mean track error (nm) for the 1975-76 dependent 
test (0), the 1977-78 independent test (4), the 


1975-78 dependent test (@), and the 1977-78 
dependent test (@). The 1977-78 dependent tes: 
(a) was the only experiment which included back- 
ward integration. 
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immeradadjusted track of Typhoon Rita in Fig. 21 @xWibits a 
Similar behavior in another westward tracking storm, but 
does very well in the longmeudindal poshtroning of She steorn. 
Adjustment of the track of Typhoon Babe, shown in Fig. 22, 
has improved the velocity forecast as well as has indicated 
Some recurvature which the TCM missed (see also Fig. 17). 
fyengon Phyllis is illustrated in Fig, 23. "Repro ve @és 
an excellent forecast of Pecunia alee not achieved by the 
unmodified TCM forecast or by the regression equations based 
only on forward integration of the TCM. The velocities in- 
dicate only a slight tendency to overcompensate for the 
Slower TCM forecast. Recurvature of Typhoon Gilda, depicted 
mimerig., 24, 1S correctly predicted, but the velocities are 
excessive. Adjustment of the track of Typhoon Wendy as in 
Fig. 25 showed slight improvement over the unmodified TCM 
positions. The interesting feature in the regression-adjusted 
path is the indication that the storm would reverse its 
direction. Within 12 hours of the indicated time, Typhoon 
Wendy did change direction in this manner, but with somewhat 


smaller displacements. 
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VI. CONGIIUS IGN S 


A three-layer, primitive-equation model (Hodur and Burk, 
1978) with one-way interactive boundaries is being tested 
at the Naval Environmental Prediction Research Facility 
(NEPRF). The objective of this research was to generate 
Statistical regression equations to adjust the TCM-predicted 
MedcKse towards the best tracks. This approach is based on 
the assumption that it is possible to adjust for systematic 
model and data errors. Development consisted of deriving 
three sets of regression equations, with two sets based only 
on forward integration of the TCM, while the last set con- 
tained predictors based on both forward and backward inte- 
gration. All TCM runs were based on operationally analyzed 
data from FNWC. Time-dependent boundary conditions provided 
to the TCM were derived from analyzed (rather than predicted) 
fields, as was the case in Hodur and Burk (1978). 

The most notable improvements occured with the regression 
equations containing predictors based on backward integration 
of the TCM. The equations with predictors having both for- 
ward- and backward-integrated positions explained the greatest 
amount of variance of any set of regression equations. Ad- 
feistments to the TCM tracks at l2-hour intervals resulted 
in predictions in the 1977-78 dependent test that nearly 
met the 7th Fleet forecast error goals. The selection of the 


Vemeeity along the y-axis in the interval -lé2 to 90 hours as 
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the predictor explaining the most variance suggests strongly 
that systematic errors at or near initialization time were 
used advantageously to adjust TCM tyack Forecasts. lt 1s 
noteworthy to observe that TCM forecasts with all data samples 
used in these experiments incurred systematic errors (see 

mig. 26). 

Tt should be noted that the weak link in the regression 
adjustment scheme is the small sample size for derivation of 
the regression equations. With the exception of the 19/5-/6 
equations, the regression coefficients have not been tested 
with an independent sample of storms. Application of the 
statistical equations in the future may not produce comparable 
results due to the relatively small number of cases used to 
derive the equations. Likewise some years have persistently 
agomalous storm tracks, as in the first sample (1975-76) 
treated here. It is expected that some increase in the 
Stability of the regression equations could be achieved by 
increasing the- number of TCM forecast cases used to derive 
the equations. This is especially true for the backward 
integration set. 

In all of the experiments, the unmodified ICM forecasts 
were difficult to beat in the case of westward propagating 
storms. However, occasional improvement in velocity predic- 
tion was noted. The typhoons used in these tests had storm 
tracks in several general categories, such as westward, south- 
east-to-northwestward, northward, or recurving. Since the 
regression equations were derived from storm cases involving 


Various combinations of these track types, it was felt that 
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these widely varying tracks and subsequent errors contri- 
buted to failures of the regression equations to make proper 
track adjustments in some cases. A possible solution might 
be to vastly increase the TCM case sample size and then 
Subdivide the cases into "regression analogs" according to 
the direction of typhoon propagation. The next step would 
be to derive a set of regression equations based on each 
Subset of storm tracks. It could be expected that each sub- 
Set of equations would have inherent adjustment biases 
towards the type of track from which they were derived. Any 
or all of these subsets of equations could be applied to each 
unmodified TCM forecast. 

On the basis of the sample examined here, it appears that 
adjustment for recurvature was handled most effectively by 
the regression equations which included backward-integrated 
tracks. It is suggested that regression equations based on 
an expanded sample of TCM cases including backward integration 
should be derived to produce better adjustments and statis- 


tical equations with greater stability. 
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APPENDEPX A: “THE TROPICAL €YCEOME MIDDIEL 


meme fae FORECAST MODEL 


The coarse-mesh version of the primitive equation model 


feveloped by Elsberry and Harrison (1971), Harrison (1973), 


Ley and Elsberry (1976), and Hodur and Burk (1978) was used 


as the basis for these experiments. The model 


mesh (2° y three-layer, dry model with one-way 


1s a coarse 


interactive 


boundaries on the north and south and cyclic east-west 


boundaries. The equations used in the model are: 


du _ 36 
ot a -L(u) + fy - M x 
av 30 
at = -L(v) - fu - M Dy 
00 _ Q 
ane =a -L(6) + TC 
JW _ 2a u e) 
OD no -M las (a) 7 Oy (7) ] 
39 = 6 C oT 
dp Pp ofp 
a one 5/16 
T= (7999) 
9?1900 
ae Oooo 7 oon erode: 
where 
0 als ee oan 
Log Sy eee Lae eye? + aye ; (wS) 





7 


(A-6) 





L(S) represents the flux divergence of any scalar quantity S. 
The meteorological symbols used above can be found in the 
mast of Symbols". 

melt ficient condition for the linear compuiwa tional —ta- 
bility of the solution for two-dimensional equations govern- 


ing simple wave motion is: 


are -707— (A-7) 
where 
C = phase velocity of the fastest gravity wave 
At = time increment 
AX = horizontal grid increment 


Computational stability requires a maximum time step of 450 
seconds for this model. Ley (1975) achieved an increase to 

890 seconds by time averaging the pressure gradient term of 

the momentum equations. Shewchuk (1977) used a 600 second 

time step for testing his 1975 cases as the relocatable grid 
was extended northward where Ax is reduced. In the TCM ver- 
Sion (Hodur and Burk, 1978) used for these experiments, the 
time step was 690 seconds unless the northern boundary exceeded 
40°N, in which case the time step was 450 seconds. 

The initial step was forward in time, while the leapfrog 
time differencing scheme was used in all Subsequent iterations. 
Friction was neglected, and the consequent storm motion was 
primarily the result of advective processes and heating. A 
Bessel interpolator was used to locate the minimum wind at 
1000 MB. Then latent heating was simulated by adding heat to 


a horizontal 7 x 7 grid centered on the minimum wind, which 


i2 





was not necessarily on a grid point. Weighting of the heat- 
ing function in the vertical was 0.3, 1.0, and 0.3 at 850, 
Boo. and 250 MB respectively. In the horizontal, the effect 
of heating is smoothed out by a Cressman weighting function 


(Haltiner, 1971) and results in a less erratic storm track. 


B. THE GRID 

The forecasts are carried out on a uniform, coarse-mesh 
ie” ) Mercatow grid true at 228 Ni ore hor zomma | eg d. aie - 
val was 205.8 km. Tne domain consisted of 32 points east- 
west and 24 points north-south. The grid was oriented so 
that each storm was initially located southeast of the center 
Omeene Grid. the vertical distribution of variables 15 shown 
in Fig. A-l. Although the variables are staggered in the 
vertical, staggering on the horizontal grid was not used in 


these experiments. The 850 mb winds were used to compute the 


@avyective terms at LOOO mb. 


C. BOUNDARY CONDITIONS 

Boundary conditions on the walls were after Hodur and 
Burk (1978). In the one-way interactive (9N) model, wind com- 
ponents normal to the boundaries are adjusted so there is no 
net divergence from the forecast domain. By integrating 


ay _ 
we Ms (A-8) 


around the boundary, the corresponding w values are obtained 
(Hawkins and Rosenthal, 1965). Distance along the boundary 


is represented by s, while We is the velocity component norma| 


as 


aes - 





100 w=0 (7) 
250 U,V Ord (6) 
400 UW) (5) 
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Fig. A-l. Vertical distribution of dependent variabies 
and pressure levels for the thr 
model (arter Harrison, 1973). 
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to the boundary with the positive direction being inward. 
Date from outside the forecast grid must continuously 
be incorporated into the boundaries. This is accomplished 
as described by Perkey and Kreitzberg (1976), and following 
their notation, the prediction equation for a variable X 


becomes 


XK (455) = X(i,d)tW(G 5) gee] Ate Cl-Ww(4,5)1 
ied 





(A~9) 


where the subscripts n and p indicate the new and the previous 


value of X. The subscript m refers to the model tendency, 
while 1s is the large-scale tendency. The weighting function 
used in this model was somewhat different from that given by 
Perkey and Kreitzberg (1976). 

The following weighting function W(i,j) was constructed 
$o that a minimum amount of noise was produced near the 
boundaries 

0.0 on the boundaries 

0.05 one grid row in from boundaries 
0.25 two grid rows in from boundaries 
0.45 three grid rows in from boundaries 
0.65 four grid rows in rom bommdaries 
0.85 five grid rows in from boundaries 


1,00 on al) cotter interior apoine:s 


Ls 





The forecast fields produced must be filtered near the 


boundaries. The filter used in this model is 


P. = (l*e)F. 4 Ge. 


j j + F. 1) (R210) 


deel j- 


where Fe icottwe= Tit teredidata at polmend. F, is the unfiltered 
data at point j, and a represents the smoothing parameter. 

The filter was applied every 40 Minutes to =the six rews amd 
columns nearest the boundaries with a = 0.5. This value of 

a yields a response function of 0 for a 2Ax wave. This filter 
ieea'sO applied over the entire grid to all the prognostic 
Variables every three hours. In this manner the forecasts 
include large-scale tendencies from outside the forecast 
domain. 

The vertical velocity at the upper boundary is equal to 
meeo and is calculated at levels 5, 3, and 1, Fig. A-l, 
through downward integration of Eq. (A-5). 

In an operational mode, the time dependent boundary condi- 
tions must be specified by a global forecast model. In these 
experiments, forecast fields were not available. Therefore, 
icing a "“perfect-prog" approach; analyzed fields taken every 


12 hours were used to specify the boundary values. 


Bee INITIALIZATION 
Analyzed fields obtained from FNWC are used as input data 
for the model. Initialization of the model is accomplished 


by calculating non-divergent winds from the stream function ». 
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ime relative vorticity, the forcing function for tie stream 
function, was obtained from the observed u and v components 


c= MS ie) - s(q) (A-11) 


The normal component of the wind was adjusted so there is no 
net inflow or outflow in the domain. The stream function was 


found from the expression 
eM eae (A-12) 


A direct solver (see Faulkner and Rosmond, 1976) rather than 
Successive over-relaxation is used to efficiently solve the 
moIsson equations in the initialization process. The non= 


divergent wind components were then calculated through 


Ie ee _ oy = 
Uy M dy? ee M aa (A-13) 


An appropriate balance between the mass and motion fields 
is achieved through partial differentiation of Eq. (A-1) with 
Besmect to x and Eq. (A-2) with respect to y. Addition of 
these equations and assuming that the time rate of change of 
divergence and gradients of the map scale factor can be neg- 


lected leads to a Poisson equation 








OV 3 U 
Be a2 1,3 3 y ee of 
Ue a (E(u It 3 at ees ye sb ere 
(A-14) 


wien can be solved uSing direct methods. 
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A. FORWARD INTEGRATION 
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The regression equations which follow were produced by 


SPSS usinjy only predictors derived from forward integration 
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B. FORWARD/BACKWARD INTEGRATION (1977-78 DATA) 

The regression equations which follow were produced by 
Sos Using nadie ors derived from forward and backward in- 
tegration of the TCM. The data sample consisted of 330 cases, 
or runs of theTCM in 1977 and 1978. Predictors in the equa- 


tions are listed in order of decreasing significance. 
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3. TCM Integrated to 48 HR, and -36 HR 
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